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ABSTRACT 

Context. Red giant stars approaching the end of the evolutionary phase of the Asymptotic Giant Branch (AGB) are, inter alia, charac- 
terised by (i) pulsations of the stellar interiors, and (ii) the development of dusty stellar winds. Therefore, such very evolved objects 
' cannot be adequately described with hydrostatic dust-free model atmospheres. 

Aims. By using self-consistent dynamic model atmospheres which simulate pulsation-enhanced dust-driven winds we studied in detail 
the influence of the above mentioned two effects on the spectral appearance of long period variables with carbon-rich atmospheric 
chemistry. While the pulsations lead to large-amplitude photometric variability, the dusty envelopes (resulting from the outflows which 
' contain dust particles composed of amorphous carbon) cause pronounced circumstellar reddening. 

Methods. Based on one selected dynamical model which is representative of C-type Mira variables with intermediate mass loss rates, 
i-C we calculated synthetic spectra and photometry for standard broad-band filters (Johnson-Cousins-Glass system) from the visual to the 

near-infrared. The synthetic photometry was subsequently compared with observational results. 

Results. Our modelling allows to investigate in detail the substantial effect of circumstellar dust on the resultant photometry. The 
pronounced absorption of amorphous carbon dust grains (increasing towards shorter wavelengths; <2abs la oc with leads to 

colour indices which are significantly redder than the corresponding ones based on hydrostatic dust-free models. Only if we account 
for this circumstellar reddening we get synthetic colours that are comparable to observations of evolved AGB stars. The photometric 
variations of the dynamical model were compared to observed lightcurves of the C-type Mira RU Vir which appears to be quite 
similar to the model (although the model is not a dedicated fit). We found good agreement concerning the principal behaviour of the 
BVRIJHKL lightcurves and also quantitatively fitting details (e.g. magnitude ranges, the amplitude decrease from visual to NIR, ab- 
^ . solute magnitudes). The analysed model is able to reproduce the variations of RU Vir and other Miras in (J - H) vs. (H - K) diagrams 

■ throughout the light cycle (ranges, loops). Contrasting the model photometry with observational data for a variety of galactic C-rich 
giants in such colour-colour diagrams proved that the chosen atmospheric model fits well into a sequence of objects with increasing 
mass loss rates, i.e., redder colour indices. 

. Conclusions. The comparison of our synthetic photometry with observational results provides a further indication that the applied 

■ dynamic model atmospheres represent the outer layers of pulsating and mass-losing C-rich AGB stars reasonably well. 

Key words. Stars: late-type - Stars: AGB and post-AGB - Stars: atmospheres - Stars: carbon - Stars: variables: general - (Stars:) 
circumstellar matter 

: ■>— j 1. Introduction First, it is well known that the more evolved AGB stars show 

/\ ■ pronounced photometric variability caused by pulsations of the 

H ' In the first paper of this series, Aringer et al. (12009] from now on stellar interiors (e.g., Lattanzio & Wood 120041 . Already in 1596, 
Paper I), we investigated synthetic spectra and photometry based David Fabricius reported the distinct light variation of a notice- 
on a grid of hydrostatic model atmospheres for C-rich red giant able red star in the stellar constellation of the Whale which he 
stars. A comparison with observational data demonstrated that named Stella Mirabilis, i.e., the remarkable star. Considerable 
these models are able to reproduce the emerging spectral energy progress has been achieved since then in the field of stellar 
distributions (SEDs) and the resulting fluxes in different photo- variability with o Cet nowadays representing the prototype of a 
metric bands quite well, at least for objects with effective tem- whole class of variable stars, the so-called long period variables 
peratures above « 2800 K. However (and not surprisingly), de- (LPVs). They are characterised by variations with amplitudes of 
viations are clearly visible for the very cool carbon stars as can U p to several magnitudes in the visual and on time scales of a few 
be seen for example in Fig. 14 of Paperl for the (J-K) colours. 10 to several 100 days. Historically, LPVs were subdivided into 
This is due to the fact that two significant phenomena, not cov- Mira variables, semiregular variables (SRV) and irregular vari- 
ered by the basic assumptions of the COMARCS models, be- ables (Lb) according to the regularity of the light change and the 
come important for the evolved stars on the Asymptotic Giant visual amplitude. As a consequence of the findings of the large 
Branch (AGB), which is the evolutionary stage of the objects variability surveys of the past (e.g., MACHO, OGLE) it is now 
discussed here. 
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considered being more appropriate to characterise LPVs by their 
pulsation mode, with Miras pulsating in the fundamental mode 
and the others in different overtone modes (e.g., Wood |2000 
Ita et al. 120041 Lebzelter & Wood 120051 While the amplitudes 
of Mira lightcurves reach several magnitudes in the visual (up 
to 10 mas as for example for^Cyg, Mattei 1997), they decrease 
further in the red (<3 mag in the /-band, e.g. Fig. 5 in Hughes 
& Wo od [19901 ) and infrared (0.5-1 mag , e.g. Fig. 7 of LeBertre 
1992). Bolometric light changes are of the order of l mag (e.g. 
Fig. 5 in Whitelock et al. [2000b . 

Second, the stars on the upper part of the AGB lose mass 
via strong stellar winds with characteristic terminal velocities of 
15 km s -1 and mass loss rates (MLRs) of typically 10~ 6 M Q yr~' 
(e.g., Olofsson l2004l Fig. 6 in Ramstedt et al. |20091 >. Not only are 
AGB stars thereby contributing significantly to the enrichment of 
the surrounding ISM (e.g., Ferrarotti & Gail 2006 1, but the dusty 
outflows - producing cool circumstellar envelopes - have also 
a considerable effect on the spectral appearance. The interaction 
of the dust grains in the outflowing material with the radiation 
field of the star leads to distinctive spectral features in the mid- 
to far-infrared range (e.g., Posch et al. 120021 Cami 2002] Molster 
& Waters 120031 Lebzelter et al. 120061 Cami & Blommaert l2006l 
as well as to a general circumstellar reddening due to the over- 
all absorption from the visual to the infrared (IR) of common 
circumstellar dust species (e.g., Fig. 1 of Wolfire & Cassinelli 
[19861 Fig. 1 of Groenewegen 120061 Andersen 120071 With in- 
creasing mass-loss rates M the SEDs as a whole become more 
and more redistributed from shorter to longer wavelengths (e.g., 
Fig.4 ofBe dijn fT9871 Fig. 1 of Sylvester et al. [T999l Fig. 5.3 of 
Glass 119991 van Loon [2007 ), the colour indices of such mass- 
losing red giants are considerably increased (e.g., van der Veen 
& Habing 1988 LeBertre & Winters [1998} and the stars appear 
reddened. 

Different dynamical modelling approaches have been used 
to model the photometric properties of evolved LPVs. 

Based on the Berlin models (Fleischer et al. 119921 Winters 
et al. 2000) for C-rich dust-driven winds, Winters et al. ([1994) 
presented theoretical light variations in the wavelength range of 
0.5 -25 /im. They discussed the influence of dust formation pro- 
cesses on the shapes of the lightcurves, as for example the asym- 
metries during the pulsation cycle caused by newly emerging 
dust layers or the superposed long-term variations over a few 
cycles due to multiperiodicity effects in the dust-forming layers. 
The authors also confronted their results with observed JHKLM 
lightcurves for selected mass-losing C-type Miras showing qual- 
itative agreement in the above mentioned points. Similar com- 
parisons of synthetic lightcurves with observational photometric 
data were later published by Winters et al. (11997) and for an 
M-type Mira by Jeong et al. ( 120031) . 

By using the Australia-Heidelberg models (Hofmann et al. 
1998 Tej et al. 2003 ) for pulsating atmospheres of M-type Mira 
variables, Ireland et al. (2004ai studied the effect of molecular 
absorption at different wavelengths in the NIR on the result- 
ing light variations and compared synthetic JK lightcurves with 
observational data of O-rich Miras with quite moderate mass- 
loss rates. Very similar results were presented in Ireland et al. 
(2004b), also here the authors showed that their models were 
able to reproduce selected aspects of the observed photometric 
variabilities (not being the focus of either paper, though). 

As described in detail in Nowotny et al. (2005a| Sect. 3), both 
modelling approaches have their shortcomings when it comes to 
consistent and comprehensive modelling of the outer layers of 
evolved LPVs - for example, the Berlin wind models are based 
on an insufficient description of the dust-free inner photospheric 



layers (no molecular opacities, grey radiative transfer) while the 
Australia-Heidelberg atmospheric models do not allow for dust 
formation and the resulting stellar wind. For the synthetic pho- 
tometry presented here we applied the combined atmosphere and 
wind models of Hofner et al. J2003I I which include these two 
points (non-grey molecular opacities and time-dependent treat- 
ment of dust-formation). In this way, they provide a detailed 
and consistent description of all the outer layers of an pulsating, 
mass-losing AGB star (cf. Nowotny et al. 120 101 below NHA10). 
In order to enable a direct comparison with observational data, 
we also chose a more subtle way to compute synthetic photomet- 
ric magnitudes compared to the above mentioned studies, i.e., 
we used real filter systems instead of adopting simplified rectan- 
gular bandpasses, and we employed correspondingly the proper 
zeropoint magnitudes for an absolute calibration. Thereby, we 
carried on with the studies on synthetic photometry started by 
Windsteig et al. ( 119971 1 and further pursued by Gautschy-Loidl et 
al. (Sect. 4 of Gautschy-Loidl 2001). The latter computed mag- 
nitudes and colour indices for a variety of different filter systems 
(among these also the one used here) based on C-rich model 
atmospheres. Most of these models were hydrostatic, but they 
also analysed a few proto-type non-grey dynamic model atmo- 
spheres. Significant improvements concerning the spectral syn- 
thesis (e.g., atomic lines) have been implemented since then (cf. 
Aringeret al. 120091 . 

The aims for the work presented here were the following: (i) 
to point out the limited applicability of the static models from 
the grid presented in PapeiT when dealing with very evolved red 
giant stars for which the assumptions of a hydrostatic configu- 
ration and no circumstellar dust are not necessarily valid any- 
more, (ii) to investigate the effects of pulsation and mass loss 
on the spectral and photometric appearance of AGB stars with 
the help of one selected state-of-the-art dynamical model atmo- 
sphere (Sect.[2]l, and (iii) to test the quality of the applied dy- 
namical models in representing the outer layers of mass-losing 
C-rich LPVs in a further and different way (compared to, e.g., the 
past spectroscopic studies; Gautschy-Loidl et al. 2004), namely 
via comparing of photometric modelling results with observa- 
tional data (Sect.0). 

For this purpose, we start with an investigation of the in- 
fluence of different opacity sources on spectra and photometry 
in Sect. 14. 11 continue with the photometric temporal variations 
(lightcurves, changes in colour indices) in Sect. 14.21 and present 
a detailed comparison with observed colour-colour diagrams in 
Sect.[43] 

2. Model atmospheres and spectral synthesis 

2.1. Dynamic model atmospheres 

The synthetic photometry presented below is based on an at- 
mospheric model which simulates the scenario of a pulsation- 
enhanced dust-driven wind as described in Hofner et al. ([2003 ). 
The specific dynamical model used (model S) resembles a typ- 
ical C-type Mira losing mass at an intermediate rate. In TableQ] 
the corresponding parameters are listed. This model proved to 
reproduce various observational results quite reasonably (e.g., 
NHA10), for details we refer to App. lA.2l (onlv available in the 
electronic edition). 

Apart from the hydrostatic initial model (which is the start- 
ing point for the dynamic calculation), we studied atmospheric 
structures of model S from three different and time-separated 
pulsation cycles - namely 0boi = [0,1], [7,8\, [14,15] - with 23 
individual dynamic phases for each cycle (cf. Fig. 8 in Nowotny 
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Table 1. Characteristics of the dynamic atmospheric model for 
a pulsating and mass-losing C-rich AGB star that was used for 
computing synthetic photometry. 





Model: 




S 


(i) 




[L 3 ] 


10000 




M, 


[M s ] 


1.0 




T„ 


[K] 


2600 




[Fe/H] 


[dex] 


0.0 




C/O ratio 


by number 


1.4 


(ii) 


R* 


[Re] 


493 






[AU] 


2.29 




log (g* [cms 2 ]) 


-0.94 


(iii) 


P 


[d] 


490 






[kms- 1 ] 


4 




k 




2.0 




Am bo[ 


[mag] 


0.86 


(iv) 


M 


[M yr-'] 


2.0- 6.5 x 10-" 




it 


[km s- 1 ] 


13-15.5 




(fc) 




0.28 



Notes. Listed are (i) parameters of the hydrostatic initial model, (ii) 
quantities derivable from these parameters, (iii) attributes of the inner 
boundary (piston) used to simulate the pulsating stellar interior as well 
as the resulting bolometric amplitude Amy, and (iv) properties of the 
resulting wind. The notation was adopted from previous papers (DMA3, 
NHA10): P, Au p - period and velocity amplitude of the piston at the 
inner boundary; f L - free parameter to adjust the luminosity amplitude 
at the inner boundary; M, u - mass loss rate and outflow velocity at the 
outer boundary; (f c ) ~ mean degree of condensation of carbon into dust 
at the outer boundary. The radial coordinates in Figs. lA.ll and lA.2l are 
plotted in units of the corresponding stellar radius R+ of the hydrostatic 
initial model, calculated from the luminosity L+ and temperature T+ via 
the relation L+ = ^nR\ ctT\. 



et al. 2005b). This allowed to study the variability throughout a 
light cycle as well as any cycle-to-cycle variations. 

Following the convention of previous papers (e.g., NHA10), 
the modelling results will be marked throughout this work by 
bolometric phases 0b o i within the light cycle in luminosity (cf. 
Fig. IA.2b ) with numbers written in italics, in contrast to the vi- 
sual phases </> v , which are usually used to designate observational 
results (with </> v =0 corresponding to phases of maximum light 
in the visual). The difficulties emerging when the two kinds of 
phase informations are to be related were discussed in detail in 
App.AofNHAlO. 

2.2. Synthetic spectra 

The outcome of the dynamic atmospheric modelling are snap- 
shots of the radial atmospheric structure during the temporal 
evolution as shown in Fig. IA.2l These served as input for the 
subsequent detailed a posteriori radiative transfer calculations. 

Based on the given r-p-structures, opacities were then cal- 
culated with the COMA code which was already described in detail 
in Paper I. Additional information concerning this opacity gener- 
ation code and the various improvements in the past years can be 
found in Aring er (US D, Gautschy-Loidl (I200D Gorfer (120051 
Crista llo et al. ( 120071 1. Leder er & A ringer (12009) . Aringer et al. 
(2009), and Nowotny et al. (2010). In accordance with the dy- 
namical model (cf. Table[TJ, we also assumed - except for car- 
bon (see below) - solar-like elemental abundances for the spec- 
tral synthesis. Aringer (2005 1 pointed out, that it is important to 
be consistent in this point. We adopted the values from Anders 



& Grevesse (1989 ), except for C, N and O where we took the 
data from Grevesse & Sauval (119941 . This agrees with our pre- 
vious work (e.g., Aringer et al. 119991 Paper I) and results in a Z G 
of =s 0.02. Starting from the adopted solar abundance pattern, the 
carbon abundance was subsequently increased according to the 
C/O ratio of the model (without changing the rest of the com- 
position). Abundances of various important neutral and ionised 
species (atomic, molecular) were then computed for all layers 
of the model atmosphere assuming chemical equilibrium for 
molecule formation and equilibrium for ionisation (see Lederer 
& Aringer 2009 for a detailed discussion and an extensive list of 
references). In addition, the depletion of the element C from the 
gas phase due to consumption by formation of amC dust grains 
was accounted for. In the next step, the resulting abundances - 
i.e. partial pressures of the neutral atoms, ions, and molecules as 
a function of atmospheric depth (cf. Fig. 2 of Loidl et al. 1999 or 
Fig. 5 of NHA10) - can be used to calculate absorption coeffi- 
cients for every point of a chosen grid in wavelength A and under 
the assumption of LTE. All necessary opacity sources were con- 
sidered: 

- continuous, 

- lines (= atoms + molecules), and 

- dust. 

For the continuous absorption all the sources listed in 
Lederer & Aringer (2009 1 were taken into account. We also fol- 
lowed their approach concerning the atomic lines with line data 
taken from the VALD database (Kupka et al. 2000). However, for 
the sake of computing time only Doppler profiles were adopted 
for the description of the line shapes (see Sect. 2.4 in NHA10 
for details) instead of full Voigt profiles. Test calculations have 
shown that this simplification has only a marginal effect on the 
synthetic photometry (e.g., differences in the second decimale 
place for the V-magnitude) while the computing time is smaller 
by a factor of two. For the microturbulence velocity we assumed 
£ = 2.5 km s -1 , which is in line with our previous work (Aringer 
et al. 119971 1 19991 1 and also with observational results for AGB 
stars (e.g., Smith & Lambert 119901 Lebzelter et al. 120081 . as 
well as it is consistent with the values that were used for prepar- 
ing the gas opacity tables for the hydrodynamic modelling. For 
the molecular contributions we did not use line listsS as in 
Paper I, but used the same opacity tables^ as Gautschy-Loidl et 
al. ( 120041 ). accounting for the following C-bearing molecules: 
CO (G oorvitch & Chacke rianll994) . CN (J0rgensen & L arsson 
[19901 ), C 2 (Querci et al. [19741 ), CH (J0rgense n et a l. [19961 
HCN (J0rgensen et al [19851 C 2 H 2 (Joreensen IT99Tb . and C 3 
(J0rgensen et al. I19891 ). An overview of the individual spectral 
features of the different molecular species can be found in Fig. 4 
of Paper I or in Gautschy-Loidl et al. (2004 Fig. 1). The sig- 
nificant contributions in the diverse wavelength regions are also 
identified in the upper panel of Fig.Q] Finally, the absorption due 
to dust grains of amorphous carbon was - in consistency with the 
preceding modelling of the atmospheric structures - calculated 
using the data given in Rouleau & Martin dl99 It set AC, cf. bot- 



1 The set of molecular species included in the current version of the 
COMA code and the corresponding references of the line lists can be 
found in Lederer & Aringer (120091 . with recent updates described in 
Aringer et al. ( 12009) . 

2 Mainly from the SCAN data base (J0rgensen 1997] see also 
http://www.astro.ku.dk/~uffegj/), a detailed description of which can 
be found in J0rgensen et al. (1993 2001}. For a discussion of the 
original OS data and revisions we refer to Gautschy-Loidl (120011 
Sects. 2.1+2.2). 
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torn panel of Fig.[3jl and assuming the small particle limit of the 
Mie theory to be validQ 

Summing up the opacity contributions of all sources, we 
computed combined opacity tables k v (A,t) in the wavelength 
range of 0.2-25/iin with a spectral resolution of R- 18 000 
(as determined by the above listed molecular opacity tables). 
On the basis of these the radiative transfer (RT) was solved for 
every wavelength point A, in spherical geometry (e.g., Yorke 
1980 No rdlundfl 9841 > wit h the h elp of our dedicated RT code 
(Balluch 1986 Windsteig 1995), which was originally devel- 
oped for the studies presented in Windsteig et al. (119971 ). Note 
that no velocity effects were taken into account in the a posteri- 
ori RT (Doppler shifts), in contrast to the recent work published 
in Nowotny et al. (1201 01) . 

The resulting spectra served then as input for calculating 
photometric magnitudes as described in the following Sect. 12. 31 
With the spectral resolution of R = 18 000 there are statistically 
meaningful numbers of points (a few x 10 3 ) in the wavelength 
ranges of the chosen broad-band filters. However, the resolu- 
tion of the synthetic spectra is much lower than what would be 
needed for an appropriate description of individual lines (the typ- 
ical line widths require a sampling with R larger than ^100.000). 
Thus, our approach works similar to the opacity sampling (OS) 
methodQ Like OS spectra, our computed spectra have a statis- 
tical nature and need to be averaged over a certain number of 
wavelength points 20-100) for further spectroscopic investi- 
gations (e.g., comparison with observational data). The influence 
of such a decrease in spectral resolution on absorption features is 
demonstrated in Fig.Q] its severeness changes with wavelength 
(compare, e.g., in the middle panel A = 2.4 /mi and 3.3 /mi). In 
this case, realistic representations of observed spectra have reso- 
lutions below R= 1000 (yellow, black). For the plots contain- 
ing spectra below (Figs.[3j}^6j, the original spectra were aver- 
aged over 50 points leading to a reduced spectral resolution of 
R - 360. In addition, the sampling rate of the output wavelength 
grid should be higher (roughly by a factor of « 5 -10) than the 
nominal resolution of the spectrum with decreased resolution, in 
order to obtain a smooth spectrum (similar to spectrographs ; cf. 
sampling theorem). In the bottom panel of Fig.Q~| the effect of 
this oversampling on the appearance is illustrated by a compar- 
ison Of a low -resolution Spectrum With ^nominal = ^sampling = 900 

(black) with the corresponding spectrum having /? n ominai = 900 
and Sampling = 9000 (magenta). 



Influence of spectral resolution on molecular features 



3 The influence of using different dust data (opacity data sets, grain 
densities, etc.) when modelling C-type LPVs with dusty outflows was 
already discussed in Andersen et al. 2003 (. There the effect on the re- 
sulting winds (mass loss rates, wind velocities, condensation degrees) 
as well as on the resulting colour indices was demonstrated. 

4 In the OS tec hnique (e.g., Ekberg et al. |1986l Alexander et al. |19891 
J0rgensen 1992 ), the mono-chromatic molecular absorption is com- 
puted at a limited number of wavelength points, distributed over the 
spectral range considered (more or less randomly, but not in such num- 
bers that the spectrum is resolved). Opacities at the selected wavelength 
points are directly calculated from line lists and subsequently used to 
solve the radiative transfer. Due to the statistical nature of the OS ap- 
proximation, it is necessary to average the OS spectrum over a number 
of frequency points in order to produce realistic spectra and reproduce 
observed features (e.g., Gautschy-Loidl 2001. Aringer et al. 2009]). As 
a consequence, the OS technique is in general only applied for spectro- 
scopic studies with low spectral resolutions. 



R= 18000 
10000 
3500 
1000 
360 




"III I' 111 II' III llllll'il 

R= 18000 




0.750 0.755 0.760 

A [/im] 



0.765 



Fig. 1. Synthetic spectrum based on the hydrostatic and dust- 
free initial model {R= 18 000, grey) calculated as described in 
Sect. 12.21 and some spectra of lower resolutions (as given in the 
legend) to illustrate the effect of the rebinning on the appearance 
of the molecular features. In the top panel, the whole spectral 
range is shown, while the middle panel show the details of a 
smaller sections. The dashed line marks a calculation where all 
line contributions were neglected (theoretical continuum level). 
The bottom panel shows the difference for a spectrum of lower 
resolution if oversampled by a factor of 10 (magenta) or not 
(black). 



2.3. Synthetic photometry 

To compute magnitudes on the basis of the synthetic spectra we 
needed filter transmissions as a function of wavelength to con- 
volve with and associated calibration standards (i.e., photometric 
zeropoints). From the large number of photometric systems in 
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use (e.g., Bessell 2005) we chose standard broad-band filters to 
facilitate the comparison of modelling results with observational 
findings. Thus, we utilised as in Paper I the response curves 
R(A) e [0, 1] for the following filters: Johnson-Cousins UBV and 
Cousins RI as given in Bessell ( 119901 ), and JHKLL'M from 
the homogenized Johnson-Glass system described in Bessell & 
Brett ( 1988 ). The data were retrieved in electronic format from 
theAsiago Database on Photometric Systems (see Moro 
& Munari 2000). Bessell et al. (1998 ) quote the corresponding 
absolute fluxes /? and zeropoint magnitudes m° for the stan- 
dard Johnson-Cousins-Glass UBVRI JHKLL'M system. Those 
were derived from models of Vega (Castelli & Kurucz 1994 1 
and Sirius (as for example in Cohen et al. 1992 1 and observed 
colours of these stars. 

The magnitude m x for a particular filter X out of the 
above mentioned set is calculated by means of (see footnote of 
Table A2 in Bessell et al. [19981 

m x = -2.5log fa x -21.1 -m° x , (1) 
whereby 

&" MA)R x (A)dA 

A* s %r. (2) 

f. x " R x (A)dA 

and fx is the energy flux per unit area and wavelength interval. 
The lower and upper cut-off wavelength of the filter passband are 
represented by the quantities Axj and A Xu , respectively. Setting 
A* = fix in Ec l- [Hl eads to m x = 0. 

The absolute fluxes refer to values measured at the Earth. 
Our synthetic spectra specify the flux at the stellar surface F*(A). 
To deduce absolute magnitudes we have to "shift" the star to 
a distance of d=10pc. Ignoring the effects of the interstellar 
extinction, we get 

m - M - 5 log d [pc] -5=0 

and, thus, the apparent magnitude m equals the absolute mag- 
nitude M. The stellar flux ft at a distance d emerges from the 
identity 

L A = AttR 2 F* = And 2 f* 

with R depicting the stellar radius. 

With the prescription given in Eq. [2] we calculate f* x from 
ft . We choose to infer the magnitude by setting f* x into relation 
to the zeropoint flux /?„ rather than inserting f* x directly into 
Eq.rjJ That way the additive constants cancel out, i.e., 

m x = -2.5 log f* x - 2 1 . 1 - m° x - (-2.5 log /j^- -2U-m£) , 

=o 

and we arrive at 

M x = m x = -2.5 log iflxlflx)- (3) 

The desired synthetic absolute magnitudes (e.g., Fig.[7j 
could then be computed with the help of Eq.[3] where absolute 
fluxes /?„ for a fictitious AO reference star were taken from the 

literature (Table A2 in Bessell et al. 1998 Q and the filter fluxes 
of our model f* x were determined by convolving the synthetic 
spectra with the corresponding filter transmission curves R X (A) 
according to Eq.|2] 

5 Note that the magnitude zeropoints zp(f^) and zp(f v ) in this table 
are erroneously reverse labelled as stated by Bessell (2005 I. 



For test purposes we computed synthetic photometry in the 
way sketched above for the grid of hydrostatic COMARCS mod- 
els presented in Paper I and compared the values with the pho- 
tometry given there which was calculated following the more so- 
phisticated approach described in Girardi et al. (2002i. Averaged 
over 738 models we found mean differences of « 0^024 -0^077 
for the filters VRIJHKL. For the M filter we found a larger dif- 
ference of 07*1 3, while no direct comparison could be carried out 
for the B filter due to the missing wavelength coverage of the 
spectra in Paper I. Note that these differences are largely caused 
by a systematic shift in the zeropoints, while the standard devi- 
ations of the magnitude difference distributions for most filters 
are below «(T009. 

A note concerning the scaling of C2 opacities and the result- 
ing shift in (H-K): As described by Loidl et al. (120011 ). the gf 
values of the C2 lines in the list of Querci et al. ( 119741 ) were 
scaled down beyond 1.15 fj.m for the preparation of their OS ta- 
bles. According to Loidl et al. this should provide a more real- 
istic description of the C2 absorption compared to the unsealed 
Querci list. In Paper I we found that the applied scaling leads to 
moderate changes of the spectra in the region of 1.3-2.1 //m, 
i.e., mainly in the H-band. From a comparison of observations 
of C stars with synthetic photometry based on hydrostatic mod- 
els and computed with scaled C2 data (see Fig. 15 in Paper I) it 
was suspected that this correction leads to worse {H-K) colour 
indices. As a rule of thumb, it was recommended in Paper I to 
add O^l to the synthetic (H-K) colours of the hydrostatic grid 
to compensate for this alleged deficient scaling (equal to using 
the original line list). This was done for the hydrostatic models 
in Figs.[l5]and[T6] For the spectral synthesis based on dynamic 
models presented here we used the same OS tables as Loidl et al. 
(2001). However, no investigations on the differences in (H-K) 
colours were made for the dynamic models as the effects of pul- 
sation and circumstellar dust (cf. Fig.fTTTi are significantly larger 
than arguable changes due to uncertain molecular line lists. 

To allow for a direct comparison with the modelling results, 
all observational data (Sect.|3]l were, if necessary, transformed to 
the Johnson-Cousins-Glass or Bessell photometric system which 
was used for the calculation of synthetic magnitudes as just de- 
scribed and should serve as the standard system throughout this 
work. 

3. Observational results for comparison 

In order to test the ability of our models to reproduce the pho- 
tometric appearance of real C-rich LPVs we used observational 
data from the literature. A valuable source in this respect was 
the compilation of near-IR data for a large sample of galac- 
tic carbon-rich variable AGB stars presented by Whitelock et 
al. (|2006). They distinguished between Miras with clearly peri- 
odic lightcurves and peak-to-peak amplitudes in K larger than 
0.4 mas (i.e., their variability type In) and other LPVs with vari- 
ations that are not periodic or have amplitudes less than 0.4 mag 
(i.e., their class 2n for non-Miras). On top of the photometric 
time-series in JHKL, they also derived for a subsample of C-rich 
Miras other quantities which are useful for a detailed comparison 
with modelling results: estimates for the distance of the targets, 
the interstellar reddening in the corresponding direction, bolo- 
metric magnitudes and corrections, or estimates for the mass- 
loss rates. 

We chose one specific object from the sample of Whitelock 
et al. (120061 ) - namely the C-type Mira RU Vir - for a more 
detailed comparison as the known properties of this star (listed 
in Table[2) are roughly comparable to the model parameters 
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Table 2. Observed and derived properties of the C-type Mira 
RU Vir collected from the literature. 



Refs. 



L 




6888 


1 




[mag] 


4.94 


1 


#0 


[mag] 


1.79 


1 


distance 


[pc] 


910 


1 


(m-M) 


[mag] 


9.80 


1 


A v 


[mag] 


0.08 


1 


spectral type 




C8,le(R3ep) 


2 


Pk 


[d] 


444 


1 


P 


[d] 


433.2 


2 


AV 


[mag] 


5.2 


2 






6.7 


3 


aa: 


[mag] 


0.98 


1 


M 


[M yr'] 


2.0- 8.3 x 10- 6 


1,4,5,6 


^exp 


[kms- 1 ] 


16.9-18.4 


4,5,6 



References. (1) Whitelock et al. 
(5) Loup et al. 0993}; (6) Knapp & Morris {1985). 



2009 >; (3) AAVSO data (V m 



(2006); (2) GCVS, Samu s et al . 
in ); (4) Bergeat & Chevalier ([20051 ; 



(TableQ} and it showed a photometric behaviour quite similar 
to our model S. However, model S was not intended to be a spe- 
cific fit for RU Vir. The difficulties arising when we aim at re- 
lating dynamic model atmospheres to certain observable targets 
were discussed in Nowotny et al. (2005b). Crucial for the choice 
of RU Vir was, in addition, the fact that it is one of the few LPVs 
for which also time-series photometry in the visual and red is 
avaible (although obtained at a different interval in JD than the 
Whitelock observations). The respective BVRI lightcurves were 
published in Eggen (|1975al ) and analysed by Eggen dl975bl ). 
Furthermore, the quite low interstellar extinction in the direction 
of RU Vir (Tabled is convenient. 

From Whitelock et al. (2006 ) we adopted for RU Vir a se- 
ries of 46 data points in JHKL over a time span of « 17 years 
(i.e., =a 14 pulsation periods). The observed magnitudes were 
not transformed as the differences between the SAAO photo- 
metric system and the Bessell filters we used for the modelling 
(Sect. l2.3l ) are almost negligible according to Bessell & Brett 
JT9881 Tabl e I)B 

Eggen d!975al ) lists 27 measurements in BV as well as 23 
in RI over an interval of »2.5 years (i.e., about 2 periods). 
The BV magnitudes are given in the Cousins system (Eggen 
I1975cl ). which is essentially identical to BV in the Johnson sys- 
tem (Bessell & Weis 1987), and could, therefore, directly be 
used here. For a comparison of the RI magnitudes with our 
synthetic photometry we needed to transform these from the 
(Eggen-)Kron system to the Cousins system with the help of the 
following conversions which are based on the data of Bessell 
& Weis ( 1987 ) and were kindly provided by M.S. Bessell (priv. 
comm.): 



Ic 
R c 



I K - 0.227 + 0.0367 (R - I) K 

R K - 0.1 13 - 0.0609 (R - I) K + 0.33 (R - 
- 0.085 (R-I) 3 K 



Our aim was to combine the available photometric obser- 
vations of RUVir - coming from different light cycles sam- 



6 Especially when compared with the dynamic effects we focus on 
in this work. While, for example, the difference between (J- K) m i t and 
lg (cf. Figll lb. the transformation would be 
-0.005. 



10 



12 



14 



? f. 



• .* - - 



0.5 

1.0r 

1.5 L 

2.0 
2.5 



(J -K) iya amounts to 2 - 3 m 

(7-X')BB = 1.0(y-A')sAAO 



2000 4000 6000 

JD-2445800 [d] 

Fig. 2. Light variations of the C-type Mira RUVir in the vi- 
sual (AAVSO data; black dots) and in the K-band (adopted from 
Whitelock et al. l200"6l red points). The arrow marks JD 2447336 
which was determined to be a maximum phase (</> v =0) from the 
AAVSO data and then adopted to estimate phases for the NIR 
data (see text). On top of the regular variation some long-term 
trend can be recognised (cf. Fig. 2b in Mattei |T997) . 

pled irregularly and insufficiently (Fig. [2]) - into one composite 
light cycle with </> v e [0,1], leading to well-sampled lightcurves 
as a function of phase (Fig.[8](. However, this kind of phase in- 
formation is not given by Whitelock et al. (2006) for the NIR 
data and the values listed by Eggen d!975al) appear doubtful 
as the light maximum clearly deviates from </> v =0 (cf. Fig. 2 in 
Eggen IT975bl ). Therefore, we used the photometric data in the 
visual available for RUVir in the AAVSO database^ (Fig.|2]or 
Fig. 2b in Mattei [19971 ) to derive visual phases ourselves. With 
the help of Period84, a software to analyse photometric time se- 
ries (Lenz & Breger 2005 ), we made a Fourier-fit to the AAVSO 
V-band lightcurve in the range where the observations of Eggen 
d!975al ) and Whitelock et al. ( 2006) were carried out (the de- 
rived period of 434 d is quite similar to the values in the liter- 
ature, cf. Table[2](. Thereby, we could determine the points in 
time of maximum phase for both time series to JD 2441262 and 
2447336, respectively. This allowed to compute phases (p Y for 
each measurement and plot the latter co-phased as shown in 
Fig. [8] Subsequently, we computed sinusodial fits (one Fourier 
component) to the resulting combined lightcurves in all filters, 
again using Period84. The fits are shown in Fig. [8] too. 

Apart from the time-series photometry we also used data 
from a few other sources. The observational data adopted from 
Bergeat et al. (2001 1 Table 4) in Fig. [15] were compiled by those 
authors from the literature. According to Knapik et al. (1997]) the 
measurements were obtained in the "Arizona system or close to 
it" (i.e., Johnson), and can be directly compared (Bessell & Brett 
1988, Table I). 

In addition, a representative sub-sample of the extensive grid 
of hydrostatic model atmospheres of Paper I (covering the whole 
range in colours) was selected and the corresponding photometry 
utilised for Fig. [TBI 

4. Modelling results 

4.1. The extended atmosphere and the role of dust 

The differences between the relatively compact hydrostatic and 
dust-free initial model and the developed model with dust par- 



7 http://www.aavso.org/data/lcg 
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tides occuring in the stellar wind (Fig. IA.2t shall be discussed 
in the following, as well as the effects of the different opac- 
ity sources on spectra and photometry. In Fig. [3] we compare 
synthetic low-resolution spectra (Sect. l2.2l i based on the initial 
model and one phase of the dynamic model for which the differ- 
ent sources of opacity were taken into account in the a posteriori 
radiative transfer or not J| Our spectra cover the wavelength range 
where the bulk of the flux is emitted by the star, with the SED 
peaking around 1-3 //m. Across the whole range significant con- 
tributions of (atomic and) molecular features are found, which is 
characteristic of late-type giants. On top of the absorptiorQ due 
to molecules, there is also the severe effect due to amorphous 
carbon dust in the outflow (> 2/?*; cf. Fig. IA.2f ) clearly visible. 
The absorption features in the visual and red (mainly CN, C2) 
appear attenuated, this can become more pronounced for other 
phases of the pulsation cycle as can be seen in Fig. [6] Even more 
noteworthy is the overall shift of the SED towards longer wave- 
lengths. As the circumstellar amC dust grains efficiently absorb 
radiation of the star in the visual (oc Q^Jy) and re-emit it in the 
IR (oc 2y bs ^v(!Td)X the flux is considerably redistributed. In this 
way, the dust is not only the trigger for the stellar wind, but has 
also a pronounced effect on the spectral appearance. 

This noticeable circumstellar reddening can also be recog- 
nised in Fig.|4] There the photometric magnitudes for the cho- 
sen standard filters (Sect. l2.3i based on the spectra of Fig.|3]are 
plotted. It is interesting to see that only small differences can 
be found between the SED based on the hydrostatic model at- 
mosphere and the SED for the dynamic phase for which only 
continuous and line contributions (atoms, molecules) were taken 
into account while dust absorption was neglected. On the other 
hand, a significant dimming of a few magnitudes in the visual 
can be seen if the dust opacity is accounted for in the a posteri- 
ori radiative transfer (consistent with the dynamical modelling). 
Increasing with decreasing wavlengths this dimming reflects the 
characteristic absorption of amC grains shown in Fig. [3] (inset 
of the bottom panel). This fundamental behaviour is valid for 
every phase as can be concluded from the bars in Fig. |4] depict- 
ing the range of variation of the synthetic photometry during the 
light cycle. While the amplitude in the V- filter only amounts to 
AV w l mag if the dust is not accounted for, we find x:4 mag for 
the full calculation including dusty effects. Note, that the ampli- 
tudes in Fig.|4]are slightly smaller than the ones plotted in Fig. [9] 
as the former only cover one pulsation cycle while for the lat- 
ter all three cycles of Fig. [7] were accounted for (cycle-to-cycle 
variations). 

The fact that the dust absorption in the optical is more rele- 
vant than the effect of the modified atmospheric structure (more 
extended, changing with time; cf. Fig. IA.21 i on the resulting 
molecular featureo - which could presumably play an impor- 
tant role as well - is illustrated in a different way in Fig. [5] If 



8 The latter is not consistent with the dynamical modelling (no con- 
servation of energy, i.e. not the same total flux) but useful for illustration 
purposes. 

9 Concerning the molecular features seen in emission around » 7fxm 
and w 14^jm we refer to the discussions in Gautschy-Loidl (2001 
Sect. 3.3) and Gautschy-Loidl et al. (2004, Sect. 6), as this effect is not 
relevant for the photometric results in this work. 

10 It may be the other way round for M-type LPVs where the cir- 
cumstellar dust is probably more transparent in the visual (e.g., Woitke 
l2^06H2007l H6fner l2008l Bladh et al. lZOTTb while the even larger pho- 
tometric variations over the light cycle in this wavelength region are 
caused by the pronounced changes in absorption due to the strongly 
temperature sensitive TiO features (e.g., Lattanzio & Wood 2004] Reid 
& Goldston l2002l Hronf2006). 



only the atomic and molecular lines are taken into account the 
dynamic model varies in this diagram around the hydrostatic 
model (not being equal at any phase, though) with deviations 
of approximately + (F2 in the colour indices. Switching on the 
dust absorption leads to a quite elongated loop which is offset 
by more than l mag in both colours. Although the influence of 
molecular features are apparent, the effects due to dust are even 
more prominent. 

4.2. Temporal variations 

For all the modelling results presented below all opacity sources 
described - i.e., atomic and molecular lines as well as dust - are 
accounted for. 

Low-resolution spectra for model S are shown in Fig.|6] il- 
lustrating the temporal variations of the spectral appearance of 
a pulsating C-type Mira with dusty outflow (see Fig.A.l in 
NHA10 for the corresponding plot of a pulsating model atmo- 
sphere without wind). Coupled to the sinusoidally varying lumi- 
nosity input at the inner boundary of the model (cf. Fig. IA.2b ). 
we find significant spectral variations during the light cycle. For 
the radial structures we find in Sect. lA.2l that the dynamic phases 
never resemble the hydrostatic atmosphere, a fact that also shows 
up here for the resulting spectra. The spectrum based on the hy- 
drostatic initial model is different from the corresponding spectra 
of the dynamic model with the fully developed wind at any phase 
(pboi. The spectra show a change in the balance of the whole SED 
(0boi —0.41 vs. 0.63) and a significantly modified appearance 
of the molecular features with a pronounced attenuation below 
2/im. 

In Fig. [7] the photometric variations in the chosen standard 
broad-band filters for the three pulsation cycles of model S we 
analysed are shown. As can be seen, the dynamic model be- 
comes fainter in comparison with the initial model for the fil- 
ters B VRIJH, while it appears brighter for the filters L'M. In the 
K-band, the result of the dynamic calculation oscillates around 
the hydrostatic magnitude. The corresponding photometric data 
for RUVir are shown in Fig. [8] where the observational data 
are plotted co-phased (cf. description in Sect. [3} and repeated 
for better illustration. A comparison of the synthetic lightcurves 
with the observed time-series photometry reveals that our model 
describes the basic behaviour of C-type MirasQ reasonably well. 
The continuous increase in brightness with longer wavelengths 
(from B to L) is reproduced; note the quite similar range in 
magnitudes on the ordinates. We also find similar absolute mag- 
nitudes compared to the model if the apparent magnitudes of 
RU Vir are shifted by the distance modulus listed in Table|2] 

Since the observed photometry is adopted from different 
light cycles and then combined, the scatter around the sinusoidal 
fits in Fig.|8]reflects cycle-to-cycle variations. For the model one 
can recognise such differences between different cycles, espe- 
cially in the visual and red, directly from the overplotted individ- 
ual lightcurves in Fig. [7] The cycle-to-cycle variations are more 
distinct towards shorter wavelengths and may be attributed to 
dust formation of varying intensity during the different cycles of 
the model. 

The observations of RUVir show the known trend of IR 
lightcurves lagging behind the visual lightcurves by « 0.1 -0.2 
in phase (see also the extensive discussion concerning this topic 
in NHA10, App. A). For the modelling results we find a phase 
shift the other way round, with light minima in the visual de- 
tectable later in phase by A^boi ~0.1. In NHA10 we found the 

1 1 assuming RU Vir to be typical 



7 



W. Nowotny et al.: Synthetic photometry for carbon-rich giants II. 



same effect based on a different dynamic model atmosphere 
(model M), and we showed that both, the amplitude and the 
shape of the visual light curve are dominated by dust opacities. 
The distinctive attenuation of molecular bands in the visual and 
red wavelength range of our current synthetic spectra for phases 
near to or shortly after the light minimum (Fig.|6) fits into this 
scenario. For model M in NHA10, the phase of the visual min- 
imum is closely linked to the emergence of a new dust layer, 
whereas the NIR light curves are close to the (sinusoidal) bolo- 
metric light curve (determined by the variable luminosity at the 
inner boundary), both in shape and phase. Therefore, the dis- 
crepancies in phase lags between the visual and NIR compared 
to observations most likely indicate problems with the dust de- 
scription in the models, which has two basic aspects: (i) the time- 
dependent formation and growth of grains, described by the so- 
called moment method, and (ii) the grain opacities, described by 
the small particle limit of Mie theory. Dynamic models indicate 
that typical grain sizes in C-rich AGB stars may exceed the range 
where the widely-used small particle limit is a good approxima- 
tion (e.g. Winters et al. 1 19971 Mattsson et al. 120101 Mattsson & 
Hofner 1201 11 1 which may have consequences for the synthetic 
spectral energy distributions and light curves. A detailed discus- 
sion of this point will be presented in a future paper. 

From observational studies of evolved LPVs it is known that 
the variability is, in general, strongest in the blue/visual (e.g., 
Eggen ll975bl Mattei ll997l Mennessier et al. 1 19971 1 and s lightly 
less pronounced in the red (e.g., Fig. 5 in Hughes & Wood 1990 
Lockwood & Wing |1971| l, while the photometric amplitudes 
decrease towards the NIR (e.g., Fig. 7 in LeBertre 1992 and 
LeBertre 1993] Percy et al. 2008). This fact is also apparent 
for RUVir in Fig. [8] The dynamic model atmosphere used for 
our modelling reproduces the characteristic decrease of ampli- 
tude with wavelength^ quite well as can be seen in Fig. [7] A 
direct comparison of the photometric amplitudes in a quantita- 
tive way can be drawn from Figsj9]and[T0] In the observational 
plot we also included the value for AV from the GCVS catalogue 
(Table|2|, as the rather limited time-series of photometric data in 
the visual obtained by Eggen (11975al l may not cover the full am- 
plitude as it can be found for studies over a large number of light 
cycles (cf. Fig. 2b in Mattei 1997). The reason for this are the 
known cycle-to-cycle variations which are also responsible for 
the Fourier amplitudes always being smaller than the total am- 
plitudes of individual measurement (»? max -m m m). The fact that 
the Fourier fit underestimates the amplitude can be recognised 
in Fig. [10] for the NIR filters. Our model also shows the corre- 
lation AK « Amboi (Fig. [9]) which was found observationally for 
Miras by Whitelock et al. (120001 their Fig. 5), although this was 
investigated only for O-rich targets with moderate MLRs. 

Figure[T2] shows the variations of RUVir in various colours 
based on the photometry shown in Fig. [8] The known trend 
(e.g., Payne-Gaposchkin & Whitney 1976, Lebzelter & Posch 
120011 ) of bluer colours around light maximum and redder colours 



12 This general qualitative trend (mainly determined by the dust ef- 
fects; cf. Fig.|6]or Fig.A. 5 in N HA10) was already reproduced by the 
models of Winters et al. (1994, 19971. Note, however, that their grey 
models for LPVs with heavy mass loss have optically thick dust shells 
and Winters et al. H99H could only compare synthetic lightcurves with 
IR data of dust-enshrouded stars which are very faint in the optical. The 
dynamic model atmospheres applied here are also suitable to describe 
Miras of intermediate mass loss where the inner dust-free photosphere 
can be observed spectroscopically (molecular features). This allows a 
comparison of the lightcurves from the visual to the IR (Fig.[7} with 
the corresponding data of objects for which the light variations in the 
optical can be observed (Fig.[8j. 



around minimum can be recognised, at least for selected colour 
indices, such as (J-K). Still, we are facing a severe scatter, re- 
sulting from the merging of measurements from very different 
periods into one combined light cycle. The variation can be seen 
much clearer for (V-I) and (V-K), where we derived simulated 
colour curves (dotted lines) from the sinusoidal lightcurve fits of 
Fig. [8] as no simultaneous measurements of the individual mag- 
nitudes are available for these colour indices. However, the am- 
plitude in these colour indices has to be regarded as a lower limit 
(see discussion above). This is illustrated in the panel for (J-K) 
of Fig. [12] where we overplotted the original data with the same 
kind of simulated variation. While the amplitude A(J - K) of the 
simultaneously obtained data amounts to 0^91, the simulated 
colour curve shows an amplitude of only 0^47 . The correspond- 
ing temporal variations in colour indices for the dynamic model 
are shown in Fig.QT] We find a good agreement with the obser- 
vations of RU Vir in the principal behaviour, i.e., the absolute 
values and the amplitudes. Comparing for example the (J-K) 
colours in Figs.QTJandQ~2]we find the model resembling RUVir 
very well. Despite some cycle-to-cycle variations, the model ap- 
pears clearly bluer around maximum phases 0boi ~ 0.0 and red- 
der at minimum phases 0boi ~ 0.5. 

Noteworthy in Fig.QT] is the difference between the colour 
indices of the hydrostatic initial model also given there [e.g., 
(J-K) - 1.46] and the colour range covered by the dynamical 
model during a pulsation cycle [e.g., (J-K) » 2.5-3.5]. For the 
latter we find significantly larger values due to the circumstellar 
reddening caused by the grains of amC dust in the stellar wind 
(Sect. l4.ll i. This fact provides the explanation for the discrep- 
ancies identified in Paper I. There we found distinct deviations 
in various colour indices between the hydrostatic COMARCS 
models and observed carbon stars (cf. Figs. 14, 15, 16, and 
18 in PapeiT). The very evolved and cool C-type giants ap- 
pear much redder than predicted by the COMARCS models, 
which are found to deliver realistic results only for objects above 
r e ff s» 2800 K. The samples of observed stars below this approx- 
imate border have, for example, (J-K) colours larger than 2 mag 
as can be seen in Figs. 14 and 18 of Paper I. While the hydrostatic 
models there show the discussed colour reversion effect with de- 
creasing effective temperatures which leads to (J-K) colours of 
< 1 .6 mag (similar to the hydrostatic initial model analysed here), 
the dynamical model exhibits (J-K) colours redder by at least 
j mag Thus, the circumstellar reddening results in colour indices 
much closer to observed values of stars where the assumptions 
for the COMARCS models (hydrostatic configuration, no cir- 
cumstellar dust) are not adequate anymore. In other words, the 
dynamic model atmospheres with dusty outflows - as model S, 
studied in this work - represent the only self-consistent method 
to describe (colours of) mass-losing LPVs. This is even more 
true if we proceed from the sample of quite moderate carbon 
stars of Bergeat (120011 1 used in Paper I to the sample of more 
evolved (redder, higher MLRs) C-type Miras of Whitelock et al. 
(2006) used for comparison with our modelling results below 
(e.g., Fig.Q~5]l. 

4.3. Colour-colour diagrams 

In the following we want to compare the results of our synthetic 
photometry with colour-colour diagrams (CCDs) usually investi- 
gated in observational studies. We chose one type of CCD which 
is commonly used in AGB research, namely (J-H) vs. (H-K). 
In a first step, we compare model S with a selected sample of ob- 
served C-Miras (among them RU Vir) that show similar MLRs 
and (variations of) colours to demonstrate that the dynamical 
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model provides realistic NIR colour indices. In the second step, 
we contrast the synthetic colours with the observational data of 
the whole sample of C-rich Miras from Whitelock et al. (2006 1, 
which covers a wider range in parameters (especially MLRs). 

In Fig.[JJwe plot (J-H) vs. (H-K) for the three differ- 
ent periods of model S we made use of. A significant variation 
of approximately half a magnitude in both colours during the 
light cycle can be recognised. In Fig. [14] we show photometric 
measurements for selected galactic C-type Miras adopted from 
Whitelock et al. (2006) for comparison. Several targets from 
their sample were found to populate the same region (note the 
same scaling of Figs.[T3land[T4b in this CCD, also the variations 
over time in the respective colour indices show similar ranges. 
This is especially true for the reference object RUVir. All of 
the targets shown in Fig.[14]have MLRs of * 1-3 x 10~ 6 M Q yr~' 
(Whitelock et al. 2006), which is quite similar to the correspond- 
ing value of the model (Table[TJ. Eye-catching is the fact that 
the model loops through the (J -H)-(H -/T)-plane once every 
pulsation cycle, due to the phase-dependent dust absorption and 
the different behaviours of molecular features contributing to the 
emerging spectra at the respective wavelengths within different 
filters. In the model (see the loops with different plot symbols 
and colours), we also found some cycle-to-cycle variations, es- 
pecially for phases around light maximum. This reflects the fact, 
that atmospheric structure and dust formation can differ severly 
for the same phases 0boi of different cycles of the model. 

It is not possible to investigate the corresponding variations 
as a function of phase for RU Vir directly as the individual data 
points in Fig. [14] were obtained over several periods of the star 
with an insufficient sampling rate (Sect.[3] Fig.[2]l. However, we 
derived simulated variations based on the Fourier-fits of the 
lightcurves (Fig. [8}, just as it was done for Fig. [12] This re- 
sults in a loop - also shown in Fig. [14] (dashed line) - similar 
to the model, with 0t,oi ~ 0.0 at the lower left of the ellipse and 
0boi ~ 0.5 at the upper right end. The range covered by the loop 
is again smaller than the corresponding range of the individual 
measurements as the fits tend to underestimate the variations 
(cf. Sect. l4.2t . We noticed that while the model loops counter- 
clockwise, the average variation of RU Vir follows a clockwise 
movement^ Hints for the latter were also found by other ob- 
servational studies in the past (Fig. 13 in Smak [T9641 Payne- 
Gaposchkin & Whitney [19761 Fig. 2 in Alvarez & Plez IT9981 . 
A similar disagreement between dynamic models of M-type ob- 
jects (with a parameterised description of dust formation) and 
real stars was also noticed by Lebzelter et al. (120101 their Fig. 20) 
in the context of equivalent width variations of spectral lines over 
the light cycle. A more detailed investigation of this effect in the 
future may reveal interesting insights into the adequateness of 
certain modelling assumptions (e.g. the inner boundary condi- 
tion simulating the stellar pulsation). 

To put the selected dynamic model into a broader context, we 
compiled further observational data for the same kind of CCD. 
This allowed to investigate the position of model S (as well 
as RUVir) on an empirical MLR-sequence. FigurefTBl shows a 
compilation of literature data for a wider range in both of the 
colour indices compared to Fig. [14] Note that while in the lat- 
ter individual measurements at different epochs are shown, the 
colours plotted now in Fig. [15] represent mean values and every 
point stands for one object. For the targets from Whitelock et al. 
(2006) we used the data prior to dereddening, enabling us to in- 



13 This is related to the fact that the phase shifts between lightcurves 
at different wavelengths show opposite trends in the model compared to 
observations (Sect. l4.2V See also App.lBl 



elude also stars with small photometric variations (non-Miras). 
For illustration purposes there are also measurements for carbon 
stars from Bergeat et al. (20(03 as well as synthetic photome- 
try based on a number of hydrostatic model atmospheres from 
Paper I shown. Despite the minor inconsistency that the latter 
colours are not subject to interstellar reddening, the comparison 
in Fig.[l5]is useful to demonstrate the influence of pulsation and 
dust effects on colour indices in general. Late-type giants with 
C-rich atmospheric chemistry can be found over quite a range in 
NIR colours. However, the hydrostatic models (at least the rele- 
vant sub-sample with temperatures above 2800 K) are only able 
to describe the observations on the very blue end of the sequence 
where we find objects exhibiting rather moderate photometric 
variations (or none at all) and low mass loss rates (or none at all). 
Both of these properties become more and more pronounced if 
one proceeds to the upper right corner of this CCD. There we 
find distinct Mira variables with much redder colours caused by 
the significant mass loss of this objects (cf. Fig. 16 of Whitelock 
etal.[2006]l. 

This fact is demonstrated in Fig. [16] in a more quantitative 
way by using only the Mira subsample of Whitelock et al. 
(2006), for which these authors derived and listed useful quanti- 
ties (esp., dereddened colours, amplitudes of photometric varia- 
tions, MLRs). With its intermediate mass-loss rate and an pho- 
tometric amplitude of AK = 0^98 (Table©, RUVir is located 
somewhat in the middle of the sequence of increasing MLRs and 
AK which can be recognised in Fig. [16] The same applies also 
for model S (with an (M) of «4.3 x lO _6 M yr _1 ), the average 
synthetic colour of which falls nicely onto the observed sequence 
for C-type Miras. Also the amplitude of model S (AK = l mag ) fits 
quite well to what is observed for targets populating the same re- 
gion in the CCD of Fig. [16] This basic finding is independent of 
whether the interstellar dereddening was accounted for (Fig. [TBI 
or not (FigfTBTl. The difference between the mean colour of the 
dynamical model and the corresponding value based on the ini- 
tial model illustrates again the effect of circumstellar redden- 
ing. For an easier comparison of the full sequence of observed 
Miras with the temporal variations of model S in Fig. Qj] we also 
marked the smaller range of this plot in Fig. [TBI 

5. Summary 

The hydrostatic and dust-free COMARCS model atmospheres 
presented in PapeiT are successful in reproducing the photomet- 
ric properties of carbon-rich giants with effective temperatures 
down to approximately 2800 K. However, this type of models is 
only of limited use for the interpretation of observations of more 
evolved C-rich AGB stars. For such objects dynamic effects - 
namely the pulsating stellar interior as well as the development 
of a stellar wind leading to a circumstellar dusty envelope - be- 
come significant. To cover these we need to apply more sophis- 
ticated atmospheric models, including atmospheric dynamics. 

In this work we studied the influence of pulsation and cir- 
cumstellar dust on the photometric appearance of late-type gi- 
ants with the help of one selected dynamic model atmosphere 
(model S) representing a typical C-type Mira losing mass at an 
intermediate rate. For such models the wind is driven by radi- 
ation pressure on amorphous carbon dust particles forming in 
the cool outer atmospheric layers. For various points in time we 
carried out detailed RT calculations based on the correspond- 
ing atmospheric structures given by the dynamic model. The 
resulting synthetic spectra were convolved with filter transmis- 
sion curves of standard broad-band filters in the visual and red 
(Johnson-Cousins UBVRI) as well as in the NIR (Johnson-Glass 
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JHKL'M). By applying adequate zeropoints we obtained syn- 
thetic photometry for the model as function of phase during sev- 
eral pulsation cycles. 

The modelling of photometric properties can be used to 
study effects which are not accessible to observations. We inves- 
tigated, for example, the influence of different opacity sources 
on the SED. Thereby, we found that the effect of the much more 
extended atmosphere of the fully developed model (wind) on the 
prominent molecular features is comparatively low (still appar- 
ent, though) when we look at SEDs based on the photometry in 
the complete set of chosen filters. On the other hand, the char- 
acteristic absorption of the amC grains gives rise to pronounced 
circumstellar reddening and a significant change of the SED. The 
outer layers of mass-losing LPVs are qualitatively different from 
hydrostatic and dust-free stars because of the dynamic processes 
and the dust formation taking place. This fact is reflected clearly 
in the radial structures, the resulting spectra, and the global SEDs 
as described by the photometry. Once the models are success- 
fully tested against observations, they are also useful for predict- 
ing light variations in filters where no observed data are available 
or for contributing to population synthesis studies. 

Contrasting the synthetic photometry with observational re- 
sults from the literature provided important and crucial infor- 
mation. Especially, the C-type Mira variable RU Vir turned out 
to be a useful reference object for several reasons (photometric 
data in several filters available, similar to the dynamical model, 
etc.). Comparing the observed BVRIJHKL lightcurves of this 
target (Fig. [8) with the corresponding photometric variations of 
the model (Fig. [7]) revealed good agreement concerning several 
aspects: 

- the objects become, in general, brighter in magnitudes when 
scanning the filters with increasing central wavelengths (i.e., 
from B to L) 

- the range in magnitudes covered by this sequence, 

- absolute magnitudes when adopting literature values for the 
distance modulus of RU Vir, 

- occuring cycle-to-cycle variations, and 

- the trend of photometric amplitudes decreasing from the vi- 
sual to the IR. 

The circumstellar reddening caused by the dusty envelope 
leads to considerably larger values for all kinds of colour in- 
dices when compared to the corresponding data based on hydro- 
static (COMARCS or initial) models. Therefore, the dynamic 
model atmospheres are more appropriate to describe evolved 
AGB stars. The applied model S is able to reproduce absolute 
values and variations of diverse colour indices of RU Vir quite 
well. 

We also tested the model against observations by looking at 
NIR two-colour diagrams. We found several observed Carbon- 
Miras showing variations in the (J -H)-(H -/Q-plane which are 
very similar to the synthetic photometry in terms of approximate 
regions and range of variations during the light cycle. Both, the 
model and RUVir describe loops in such CCDs, although the 
orientation turned out to be contrary. Average colour indices in 
the same diagram compiled for various observed C stars show 
the known trend of redder colours with increasing mass loss 
rates. Losing mass at an intermediate rate, RU Vir can be found 
roughly in the middle of the linear sequence crossing the (J-H) 
vs. (H-K) diagram. Model S has almost the same mean colours 
and, thus, fits very well into this global scenario. 

Notwithstanding these successes of our modelling, some dis- 
crepancies become apparent from a more detailed comparison of 
the synthetic photometry with observational data (e.g. shapes of 



lightcurves, phase lags between visual and IR lightcurves, senses 
of rotation of loops described in CCDs). These are certainly 
shortcomings of the current models. However, such differences 
emerging from a strictly consistent modelling (as carried out in 
this work) can also provide unique insights into the physics tak- 
ing place in the outer layers of evolved red giants which will help 
to improve the modelling approach. 

The dynamic model atmospheres applied in this paper were 
already successful in reproducing observed mass loss proper- 
ties (M, Mattsson et al. 120101) . observed SEDs and low- 
resolution spectra (Hofner et al. 2003. Gautschy-Loidl et al. 
2004), line profile variations in observed high-resolution spectra 
(Nowotny et al. 2005a 2005bl 120101) . variations in line equiva- 
lent widths over the light cycle (Lebzelter et al. 120101 ) and in- 
terferometric measurements (Sacuto et al. 1201 It . The synthetic 
photometry presented here provides a further clue that the mod- 
els resemble the global atmospheric structures (and their tempo- 
ral variations) of mass-losing C-type LPVs to a rather reasonable 
degree. 

In the future, we plan to extend our efforts to model photo- 
metric properties and variations - carried out exemplarily and in 
detail for one selected atmospheric model, here - to the grid of 
dynamical models published by Mattsson et al. (2010). Covering 
a wide range of parameters, this will enable us to test the influ- 
ence of different parameters (e.g., MLR) on the resulting syn- 
thetic photometry. 
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Influence of omorphous corbon dust on the SED 
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Fig. 3. Illustration of the influence of different opacity sources 
(for details see text) - especially the grains of amorphous car- 
bon dust - on the resulting spectra demonstrated with the help 
of the dust-free initial model and one phase of the dynamic cal- 
culation ((f>bo\ = 0.0) of model S. Different sources were included 
in the a posteriori RT as denoted in the legend. In the middle 
panel the molecular species responsible for the major features 
are indicated. In the bottom panel (showing the same as the mid- 
dle panel but scaled linearly) the wavelength ranges of the used 
broad-band filters (« FWHM of the responses given in Bessell 
1990 and Bessell & Brett |1988j l are marked for orientation pur- 
poses, while the insert shows the absorption data for amorphous 
carbon dust from Rouleau & Martin (1991) as applied for the 
computations. 
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Fig. 4. The influence of different opacity sources on the broad- 
band photometric magnitudes demonstrated with the help of the 
initial model and one selected dynamical phase {<ph \ —0.31) of 
model S. In addition, the maximum variations throughout the 
pulsation cycle are marked with bars. Same colour code as in 
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Fig. 5. The influence of different opacity sources on the loca- 
tion in a colour-colour diagram. Shown are the hydrostatic ini- 
tial model and several dynamical phases throughout a pulsation 
cycle of model S with opacity sources accounted for as listed in 
the legend. Same colour code as in Fig. [3] Selected phases 0boi 
are denoted for orientation. 
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Fig. 6. Synthetic low -resolution spectra (plotted on logarithmic 
and linear scale, respectively) for the whole wavelength range 
considered based on the initial model as well as different phases 
of the dynamic calculation of model S. 



Dynamic model atmosphere for C-type Mira 
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Fig. 7. Synthetic lightcurves (absolute magnitudes; 0b o i = [0,7]) 
of model S computed for various filters as labelled on the right 
hand side. The data points are repeated to emphasise the peri- 
odicity. The corresponding data from two other pulsation cycles 
(0boi = [7,8] and [14,15]) are overplotted with dotted and dash- 
dotted lines. On the left hand side, the photometric magnitudes 
based on the hydrostatic initial model are drawn in for compari- 
son with the same colour code as in Fig. [3] 
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Fig. 8. Observed light variations of RUVir (cf. Sect. [3} in dif- 
ferent filters as labelled on the right hand side. The data were 
adopted from Eggen dl975al BVRI) and Whitelock et al. d2006l 
JHKL), respectively. Measurements from different periods are 
merged into a combined light cycle, each data point is then plot- 
ted twice to highlight the periodic variations. Sinusoidal fits (one 
Fourier component, see text for details) are overplotted with dot- 
ted lines to guide the eye. An analogous plot for an M-type Mira, 
namely RRSco, can be found in Lattanzio & Wood (2004; their 
Fig. 2.44). 
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Fig. 9. Amplitudes of the synthetic lightcurves for model S as 
shown in Fig. [7] (maximum value from the three different cy- 
cles minus the corresponding minimum value) plotted against 
the central wavelengths of the different filters. The amplitude of 
the variation in mb i of model S (cf. Fig. lA.2l > is marked, too. 
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Fig. 10. Amplitudes of the observed light variations of RU Vir as 
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'Wmin) of the data sample from Whitelock et al. (2006) are plotted 
with filled circles, together with the amplitude of the Fourier fit 
given by Whitelock et al. in their Table 3 plotted with open cir- 
cles. For the BVRI filters only the first kind of amplitude is avail- 
able and plotted with asterisks. In addition, the GCVS value for 
AV is shown with a filled square. The discrepancy in V is due to 
the limited sampling of the data obtained by Eggen d!975al >. 
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Dynamic model atmosphere for C-type Mira RU Vir (C-type Mira) 
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Fig. 11. Variations of different colour indices for model S, over- 
plotted with the same linestyles as in Fig.Qare the values for dif- 
ferent pulsation cycles. On the right hand side the amplitudes in 
the colour indices are listed as well as the corresponding colours 
of the hydrostatic initial model for comparison. Same colour 
code as in Fig. [3] 
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Fig. 12. Variations of different colour indices of RU Vir. The 
dotted curves were derived from the sinusoidal fits in Fig. [8] as 
no simultaneous measurements of the individual magnitudes are 
available for (V-/) and (V-K). To get an idea of how repre- 
sentative these estimates are the same simulated colour curve 
was derived for (J-K), where photometry at the same epochs is 
available as well. On the right hand side the amplitudes of the 
respective colour indices are listed (in brackets if only derived 
from the simulated variations). 
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Fig. 13. Colour-colour diagram for model S where the values for 
the three different pulsation cycles considered are plotted. For 
orientation purposes the arrow illustrates how the model loops 
through the diagram during one cycle and selected phases are 
labelled. 
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Fig. 14. Colour-colour diagram containing observations at sev- 
eral epochs for a sample of C-type Miras adopted from 
Whitelock et al. (|2006 ) as listed in their Table 2 (no correction 
for interstellar reddening). In addition, the dashed line marks 
the average variation of RU Vir during the light cycle as derived 
from the sinusoidal fits in Fig. [8] The arrow marks the direction 
of how the star passes through the loop. 
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Fig. 15. Mean colours (no correction for interstellar reddening) 
for C stars derived from the Fourier-mean JHK magnitudes listed 
in Table 3 of Whitelock et al. (2006 ) for all Miras (i.e., their vari- 
ability type In) and all non-Miras (i.e., their class 2n). Each point 
represents the time average of temporal variations (as shown in 
Fig. [14] for a subsample of the targets included here) for a spe- 
cific object. The dashed line represents the fit given in Eq. 7 
of Whitelock et al. (120061) for the Mira sample prior to de- 
reddening. In addition, the photometric data collected by Bergeat 
et al. ( 120011 Table 4) is overplotted as well as synthetic colours 
from a sub-grid of the hydrostatic COMARCS models presented 
in Paper I for comparison. 
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Fig. 16. Mean colours corrected for interstellar reddening as 
listed in Table 6 of Whitelock et al. (120061 1 for all C-type Miras 
(i.e., their class In). The dashed line represents the fit given in 
Eq. 2 of Whitelock et al. (2006 ) for the de-reddened Mira sam- 
ple. The labels along the sequence mark approximate mass loss 
rates M [M Q yr~'] (above the fit line) as well as photometric am- 
plitudes AK [mag] (below) corresponding to given colours as 
derived by Whitelock et al. Overplotted are the synthetic colours 
of model S (initial model as well as the average of several phases 
of the various cycles in Fig. [Till for comparison. The box drawn 
with dotted lines marks the range covered in Figs.[T3ltiT4l 
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Appendix A: Details of the atmospheric models 

A.I. General remarks 

The two aspects discussed in this paper - pulsation of the stel- 
lar interior and the formation of a dusty stellar wind - have a 
considerable influence on the outer layers of an AGB star and 
hydrostatic model atmospheres are, therefore, often not an ade- 
quate approach to describe the resulting complex stratifications 
and dynamic effects. As a consequence, dynamic model atmo- 
spheres were developed to simulate mass-losing LPVs (Sect.Q]). 

The majority of the AGB stars found are of spectral type M 
(e. g., No wotny et al. l2003l Battinelli et al. I2U031 Rowe et 
al. I2005l l and have oxygen-ricrQ atmospheric chemistries with 
C/O < 1 . This is also reflected in the mineralogical species con- 
stituting the dust grains that occur in the winds. Observational 
spectroscopic studies revealed a rich mineralogy in the dusty 
outflows of M-type giants (e.g., Molster & Waters 2003] and 
references therein). However, the dust formation process in the 
atmospheres of such objects is so far not fully understood (cf. 
Hofner 2009) and the open question concerning the driving 
mechanism for the wind in the O-rich case is not solved, yet 
(Woitke I2006l t j2007l Hofner 2007]l. One possibility to drive a 
wind was sketched by Hofner (2008 ) with the help of big (x: ^m) 
Fe-free silicate grains (i.e., Forsterite Mg2Si04) and their sub- 
stantial radiative scattering cross section. 

As a star evolves during the AGB phase, a combination of 
nucleosynthesis and convection processes can lead to a dras- 
tic change of the chemical composition of the outer layers 
(e.g., Busso et al. [1999 , Herwig 2005). Having been converted 
to a carbon star (Wallerstein & Knapp 1998 ) with an atmo- 
spheric chemistry characterised by C/0>1, the spectrum of 
this star is strongly changed and shows prominent features of 
C-bearing molecules (Paperl, Fig. [3). This transformation to 
spectral type C is, in addition, highly relevant for the circumstel- 
lar dust chemistry which is believed to be simpler than the cor- 
responding one of M-type stars. Apart from the few dust species 
identifiable via their characteristic emission features, as for ex- 
ample SiC or MgS (see Molster & Waters 2003), the majority 
of the formed dust grains are composed of amorphous carbon 
(amC). Although they are not recognisable in the spectra by a 
distinctive spectral feature because of the featureless extinction 
properties (e.g., Fig.[3]or Andersen et al. 11999) . amorphous car- 
bon dust plays a crucial role for driving stellar winds via ra- 
diation pressure as outlined for example in NHA10. Moreover, 
the formation and evolution of amC grains can be treated in a 
consistent way (Gail & Sedlmavr 119881 Gauger et al. 11990) in 
numerical models. This enables us to calculate detailed and self- 
consistent models for the atmospheres of pulsating red giants 
including the developing winds. 

A.2. Dynamic model atmospheres 

For the synthetic photometry presented in this work we used the 
dynamic model atmospheres presented in Hofner et al. (|2003 ). 
These models simulate pulsation-enhanced dust-driven winds 
(caused by radiation pressure on amorphous carbon dust parti- 
cles), which is the most widely accepted scenario for mass loss 
on the AGB in the C-rich case. They are well suited to describe 
the complex behaviour of all atmospheric layers - from the in- 
ner photosphere out to the cool wind region - of pulsating AGB 



14 In fact this does not mean O-enriched, but rather the standard value 
for abundances of elements in the cosmos. 
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Fig. A.l. Movement of mass shells with time at different depths 
of model S, depicting a pulsation-enhanced dust-driven wind. 
The shown trajectories follow the evolution with time of certain 
matter elements starting from the distribution of adaptive grid 
points (higher density of points at the locations of shocks) at an 
arbitrarily chosen point in time (in this case at the end of the 
plotted sequence). Colour-coded is the degree of condensation 
of the available carbon into amorphous carbon dust grains. 

variables with intermediate to high mass-loss ratesF^lThis is ac- 
complished by a combined and self-consistent solution of hydro- 
dynamics, frequency-dependent radiative transfer and a detailed 
time-dependent treatment of dust formation and evolution (for 
details see Hofner et al. 120031 ). Additional descriptions of the 
dynamical models can be found in Gautschy-Loidl et al. (2004 ), 
Nowotny (120051 ). Nowotny et al. (12005all2l7l0l) . and Mattsson et 
al. (12007II2010I ). 

In a recent study, Mattsson et al. (120101 ) computed a grid of 
C-rich dynamic model atmospheres and investigated the result- 
ing mass loss properties (e.g., mass-loss rates M, wind terminal 
velocities dust condensation degrees in the outflows f c ) as a 
function of stellar and piston parameters. Here we made use of 
only one selected atmospheric model, the parameters of which 
are listed in TableQ] This model S is representative of a typical 
C-type Mira with intermediate mass loss and was successful in 
the past in simulating observational results as diverse as low- 
resolution spectra (Gautschy-Loidl et al. 120041 ). line profile vari- 
ations (Nowotny et al. [2005a 2005b 2010) and interferometric 
properties (Paladini et al. 120091 ). 

In Fig. IA.ll the motions of layers at different atmospheric 
depths for the chosen model are illustrated, while the corre- 
sponding atmospheric structures are shown in Fig. IA.21 The dy- 

15 The mechanism for mass loss rates below = 1(T 7 M Q yr~' (SRVs), 
where the M-P-relation (e.g., Fig. 1 of Vassiliadis & Wood 1993 Fig. 22 
of Olofsson et al. [19931 Fig. 7.25 o f Olof sson l2004l no l onger holds, is 
still dubious (Gautschy-Loidl et al. 120041 Mattsson et al. l2010t . 
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Fig. A.2. Characteristic properties of model S. Plotted in panel (a) is the bolometric lightcurve resulting from the variable inner 
boundary (piston), the circles mark instances of time for which snapshots of the atmospheric structure were stored by the radiation- 
hydrodynamics code, partly labelled with the corresponding values for phase 0boi- The dotted vertical lines mark the points in time of 
0boi = 0.0 and 1.0, respectively, while the horizontal line marks L+ of the initial model which is almost equal to (L) of the dynamical 
calculation (NHA10). The other panels show the atmospheric structures of the initial hydrostatic model (thick black line) together 
with selected phases 0t,oi of the dynamic calculation during one pulsation cycle (colour-coded in the same way as the phase labels of 
panel (a)). While panel (b) shows the classical plot of gas temperature vs. gas pressure used to characterise stellar atmospheres, the 
middle and lower panels illustrate the radial structures of gas temperatures (c), gas densities (d), gas velocities (e), and condensation 
degrees of the element carbon into amorphous carbon dust grains (f). 



namic calculation starts with a hydrostatic initial model charac- 
terised by stellar parameters as given in the first part of TableQ] 
This comparably compact and dust-free atmosphere is con- 
structed in a similar way as classical hydrostatic model atmo- 
spheres for red giants, especially the COMARCS models used 
in Paper I (for a quantitative comparison see App. IA.3t . The ef- 



fects of stellar pulsation are subsequently introduced by a vari- 
able inner boundary, representing a sinusoidally moving piston 
at the innermost radial point of the model with parameters as 
listed also in TableQ] Apart from the kinetic energy input, this 
also prescribes a changing luminosity input at the inner bound- 
ary (Fig. IA.2a ). for a detailed description we refer to Eqs. 1 and 
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2 in NHA10 as well as the corresponding explanations there. 
A shock wave triggered by the pulsation emerges during every 
pulsation cycle and propagates outwards through the atmosphere 
causing a levitation of the outer layers. In the wake of the shocks 
(post-shock regions) the physical conditions - strongly enhanced 
densities at temperatures low enough to allow for dust conden- 
sation - provide the necessary basis for the formation of amC 
grains (e.g., Sedlmavr |1994) . Radiation pressure acting upon the 
formed dust particles results in an outwards directed accelera- 
tion. Subsequent momentum transfer between the grains and sur- 
rounding gas via direct collisions (e.g., Sandin & Hofner 2004) 
leads to the development of a stellar wind. The temporally vary- 
ing radial structure of the dynamic model atmosphere (Fig. IA.2t 
deviates significantly from the hydrostatic case. Not resembling 
the intital model at any point in time, the atmospheric structure 
of the fully developed mass-losing model becomes extremely ex- 
tended in comparison, with strong local variations superposed on 
the shallow density gradient. 

A.3. Comparison of hydrostatic initial models and 
COMARCS models 

In the course of our synthetic photometry we tested how sim- 
ilar the hydrostatic intial models of the dynamical calculations 
(cf. App. IA.2l ) are compared with classical hydrostatic model at- 
mospheres like the ones used in Paper I. For the parameters of 
the initial model of model S (Table[TJ no converging solution 
could be otained with the MARCS code for numerical reasons. 
Therefore, we turned to the slightly more moderate model M 
used in Nowotny et al. (2010) for line profile modelling. The 
stellar parameters of the corresponding initial model are com- 
pared to the parameters of our closest COMARCS model in 
Table lA. ll Note that there are marginal differences as the pri- 
mary input parameters are not the same for the code to compute 
initial models and the MARCS code. The structures of the two 
atmospheric models are compared in Fig. IA.3l Apart from the 
larger range covered by the COMARCS model, the radial struc- 
tures are very similar for most of the depth points. This is also 
reflected in the resulting spectra which are shown in Fig. IA.3l 
too. To quantify the differences for the different types of hydro- 
static models we also compared the synthetic photometry in the 
chosen set of broad-band filters (Sect. 12.3V The results, as listed 
in Table lA. 21 show that the deviations are quite low. This is espe- 
cially valid in the NIR. Thus, we conclude that the initial models 
- which are constructed using fewer frequency points and no 
convection description compared to the COMARCS models of 
Paper I - are comparable to other hydrostatic atmospheric mod- 
els. They represent adequate starting points for the subsequent 
dynamical calculations, where the effects of pulsation and stel- 
lar winds lead to considerable changes in the radial structures 
(Fig.E3. 

Appendix B: Looping the colour 

A certain disagreement was found in Sect. l4.3l between the mod- 
elling results and corresponding observational data. In Fig. [13] 
one recognises that the model loops counter-clockwise through 
the colour-colour plane, while the average variation of RU Vir in 
Fig-El (dashed li ne ) passes clockwise through an ellipse in this 
CCD. It can be shown by a simple test that such a change in the 
sense of the rotation can easily be introduced by a change of sign 
in the phase shift between lightcurves in different filters, which 
may itself be linked to the dust description in the modelling con- 
text, cf. Sect.l4~2l 



Table A.l. Characteristics of the hydrostatic initial model used 
as the starting point for dynamical model M (Nowotny et al. 
[20T0b and a very similar COMARCS model. 



Model: 


initial model COMARCS 


L 




7000* 


7015 


M 


[M ] 


1.5* 


1.5* 


T e g 


[K] 


2600* 


2600* 


[Fe/H] 


[dex] 


0.0 


0.0 


C/O ratio 


by number 


1.4 


1.4 


R 




412 


413 


log (#* [cms 2 ]) 


-0.614 


-0.620* 



Notes. Those quantities representing the primary input parameters are 
marked with asterisks. 
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Fig. A.3. Comparison of the atmospheric structures of the mod- 
els in Table lA. ll (upper panel) and the resulting low -resolution 
spectra based on these (lower panel). 



To illustrate this, we used the simulated JHK lightcurves 
from Fig.|8]and applied artificial phase shifts of A(p = +0.1 such 
that the phases of light maximum in the different filters occur at 
increasing or decreasing phases <f> v . The resulting light variations 
are shown in the upper panel of Fig. IB.ll 

Based on these lightcurves, we computed colour indices 
which are shown in the lower panel of the figure. For the JHK 
variations in phase (black), the object moves along a straight line 
from bluer to redder colours in the CCD. The ranges in (J-H) 
and (H-K) are determined by the amplitude differences of the 



W. Nowotny et al.: Synthetic photometry for carbon-rich giants II., Online Material p 4 



Table A.2. Differences of photometric magnitudes calculated on 
the basis of the two atmospheric models of Table lA. ll 



Filter 


( ra COMARCS ~ m initial) 




[mag] 


B 


0.276 


V 


0.183 


R 


0.118 


I 


0.057 


J 


0.005 


H 


0.025 


K 


0.012 


V 


0.042 


M 


-0.011 



individual lightcurves. If the above mentioned artificial phase 
shift is applied to the JHK lightcurves (red, green), the star fol- 
lows a loop through the CCD which covers a larger range in 
both colour indices. In addition, Fig. IB. II demonstrates that the 
sense of rotation changes when the sign of the phase shift A(f> is 
inverted. 
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Fig. B.l. Upper panel: Simulated JHK lightcurves adopted from 
Fig-ERsinusoidal fits to the observational data there) with artifi- 
cial phase shifts of A(f> = -0. 1 / / +0. 1 (red, black, and green 
lines, respectively) imposed in all three filters. Lower panel: 
Resulting variations in an NIR colour-colour diagram, colour 
coded correspondingly. The labels mark locations occupied by 
the object at certain phases (p v , while the arrows designate the 
sense of rotation of the loops. 



